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1. 


In  an  earlier  paper  [6]  a  city  distributed  uniformly  over  a 
circle  was  considered  in  a  model  for  an  n-weapon  attack  using  the 
circular  normal  damage  function  and  aiming  error  developed  by  Hunter 
[4]  .  The  expected  number  of  casualties  was  evaluated  using  the 
Circular  Coverage  Function  p(R,r)  ([1],  [2]). 

Many  damage  curves  do  not  fit  the  single-or  two-term  circular 
normal  form  of  Hunter,  or  the  elliptical  normal  form  discussed  by  Grubbr. 
[3],  These  "cookie-cutter"  curves  are  better  fitted  by  the  damage 
function 


P(r)  * 


1 


,  0<r^k 


(1.1) 


|  Qexp(-cr8)  ,  r>k,  Q£0,c^0, 
where  r  is  the  distance  from  ground  zero. 

We  consider  a  single-weapon  attack  on  a  city  having  a  circular 
normal  distribution  of  population  about  the  origin  (see  Sherratt  [7])  ; 
the  aiming  error,  as  in  Hunter's  model  [4],  is  assumed  to  be  circular 
normal  about  the  point  of  impact  (xj,  ,yi )  ,  and  the  damage  function  of  the 
form  (1.1)  .  The  Circular  Coverage  Function,  discussed  in  [6],  will  again 
figure  in  the  evaluation  of  the  expected  number  of  casualties. 


♦This  research  was  completed  at  the  Research  Triangle  Institute. 
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Civil  Defense.  Revised  at  Southern  Methodist  University,  October  1969. 


*MTMW«S?W (U?faa^ 


2 


Section  2  evaluates  the  expected  casualties  when  the  impact 
point  is  given,  and  in  Section  3  an  expression  for  the  damage 
function,  modified  to  include  aiming  error,  is  obtained.  Section 
4  estimates  the  expected  casualties  when  the  aiming  point  (Xi,Yi) 
is  given.  To  the  writer's  knowledge,  no  previous  model  has  led  to 


such  a  result  for  the  cookie-cutter  damage  function  (1.1)  . 


'•  <n* 
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2 .  Single  weapon  blast  at  (xi,yi)  . 

Let  the  density  of  city  population  be  (as  in  Sherratt  [7]) 
P(x,y)  *  P0  exp[-(xs+ya) /2b3]  .  I 

Consider  a  weapon  blast  at  (xj .y^  ,  and  let  Tx  denote  the  expected 
casualties.  Then,  following  Fig.  1, 


where 


(xi,yj.)  Impact  point. 


(x,y) 


Fig.  1 


Ti  =  J  J  P  (x  ,y)  P  (Rj)  dxdy 

_eo 

S  *  |(x,y)  I  (x-  ,)a  +  (y-yi)®Sk3j 
5'  *  |(x,y)  |  <x,y^S  j 

ti  ■  jis  p°  exp[-  ] dxdy  +  is  -  p°q  exp[- 

s 

-  c|(x-xi)8  +  (y-yj.)  3}Jdxdy 

■  P0ii  +  poQIa » 

Ii  “  JJ  exp[-(xa+y3) /2b3]dxdy, 
o 

I8  *  '  exp[-(x3+y3) /2b3  -  c|(x-xx)3  +  (y-Yi  >  '’}]dxdy 
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Let 

dj*  =  xx2+  yia  ,dx>0 

and 

ya  =  1+2  bj  c  ,yX) 

In  order  to  evaluate  I*  and  I8 ,  we  require  the  Circular  Coverage 
Function  ([2],  [6]) 

p(R,r)  -  l/(2TT)J*Jg  exp  ^1a(xa+y3)j’dxdy 

integrated  over  the  region  R  in  which 
(x-a1) 3  +  (y-b')a  ^R8, 

a'3  +  b,a  *  ra . 

Then 

Ix  =  2rrbap{k/b,  d^b)  .  (2.4) 

We  also  require  Lemma  1  and  Corollary  1  of  [6]  . 

The  corollary  proved  there  states  (in  the  present  notation)  : 

r^,exp[-Axa+2Bx-Aya+2By]dxdy  *  nA-1exp[ (B*+B*3) »J 

.{  1-P[  (2A)  (2(B3+B,a)/A)  ]j  (2.5) 

where  A>0,  and  B, B'  are  arbitrary. 

Then,  setting  u  *  x-xx ,  v“y-yj , 

I#  =  exp[-di2/(2ba)  ]  .JJ5,exp[-(ua+v3)  Y3/  (2b*) 

-  XiU/b3  -yiv/b3]  dudv 

=  (2nb3/Ya)exp(-cd13/Y8){l-p[JtY/t>,d1/<Yb)]]'  (2.6) 

on  applying  Corollary  1  of  [6]  over  the  circle  S:  u3+va^ka , 


and  5',  and  simplifying. 


Hence,  applying  (2.4)  and  (2.5)  in  (2.3),  the  conditional  expected 


casualties,  given  impact  occurs  at  (xx ,yx) ,  are  given  by 
Tx  -  (2np0Qb3/va)exp(-cd13/y3)[l-p(kv/b,d1/(by))} 

+  2TTp0b*p(k/b,d1/b)  (2.7) 

3 .  Modified  Damage  Function. 

If  aiming  error  about  an  aiming  point  (X,Y)  is  considered, 
it  is  convenient  to  obtain  an  expression  which  modifies  the  damage 
function.  Fig.  2  shows  the  model  in  diagrammatic  form. 


(X,Y)  Aiming  point 


Fig.  2 

Let  the  distance  from  the  aiming  point  to  (x,y)  be  R,  and 
let  aa  =  1  +  2crac, 

where  aa  is  the  variance  parameter  of  the  circular  normal  distribu¬ 
tion  of  aiming  error  about  (X,Y)  .  Then  we  obtain  the  modified 
damage  function  in  the  form  of 

Theorem  1.  If  the  probability  of  an  individual  becoming  casualty 
at  a  distance  r  from  the  impact  point  of  a  weapon  is  given  by  (1.1) , 
and  if  the  impact  point  has  a  circular  normal  distribution  about 
the  aiming  point  with  variance  parameter  ca ,  then  the  probability  of 
an  individual  at  (x,y)  becoming  casualty  is  given  by 


P^  (Cas)  =pfk/a,a/a]  +  (9/Qf3)  exp(-cR  8/a8)[l-p(kQr/a,  R/ (otO) )  } 

(3.1) 

Proof:  Let  S  ^  ^  =  |(x,y)  j  (x-xj ) 8  +  y-yx8)sk8} 

where  (xx,yx)  is  the  impact  point;  thus  S  ^  is  the  set  of 

(x,y)  such  that  the  probability  of  becoming  casualty  is  1,  when 

(xx,yx)  is  specified.  Then  following  Hunter's  argument  (cf.his  Eq.(10)), 
00  00 

P,  .  (Cas)  =  f  f  P,  ,  (Cas  I  impact  in  element  dxxdyx) 

<**y>  1®  i»  {x'y} 

. Pr{ impact  i &  element  dxxdyx) dxxdyx 

=JJ^  (1/  (2TT08) )  exp£  (V  (2CTa)  ){(xi-X)  a+  (yi  -Y)  a}Jdx1dy1 

(x,y) 

+JJ  (  Qexpf-c  (x-xx) a-c  (y-yi)  a] .  exp  -(1/  (2®®)  ){(xx-X)8 
(x,y) 

+  yx-Y)8}  dxxdyi 


+  yx  -Y) 8  J Jdxx  dyx 
★  * 

=  Ix+Qis  » say . 
Again,  pat  u=x-xx  ,  v=y-yx .  Then 


(3.2) 


Ix  =  JJ<,  (2ttCT8)  -1  exp |j  (208)  _1  { (u+x-X)  3+  (v+y-Y)  8|Jdudv 
=  (2tt a8)  _1exp[-  (2C3)  -1  R3]  JJ^exp |j(2aa)  "1-[ua+va-2u  (x-X) 
-2v  (y-Y)  Jfjdudv 

=  (2170®)  "lexp[-(20®)  _1R3]  .2na®  exp  £(20*)  ~xR*j  pjk/°,R/0^j 


=  p(k/a,R/CT)  . 


Likewise 


Ig  =  (2Tra8)  _1exp(-R3/(2a3) )  JJ  exp  ||-Q'8  (u8+v3) 

:  (x-X)  u-2  (y-Y)v}/(20^J  dudv 
(l/a3)exp(-cR8/a8){-p(ka/ff,  R(Oa)"1)}- 


-2 (x-X) u-2 (y-Y)v 


•  ...k;,  -y 


on  applying  (2.5)  and  simplifying. 


The  result  follows,  from  (3.2)  .  Q.E.D. 

Aiming  error  is  now  incorporated  into  the  damage  function  by 
Theorem  1,  and  casualties  from  a  single  weapon  attack  may  now  be 
estimated . 

4 .  Single  Weapon  Attack  with  Aiming  Error. 

Denote  by  P(Cas)  the  damage  function  (3.1)  .  If  p(x,y)  is  the 
city  population  density  at  (x,y)  ,  then  the  expected  number  of  casual¬ 
ties  C  from  a  single  weapon  aimed  at  (X,Y)  is  given  by 

CC  00 

c  “  Ico  loo  p(x,y)P(Cas)dxdy 

Using  the  circular  normal  urban  density  model  (2.1) ,  let 

P  -  o/a  +  2b“c,  e>0. 

Then  Theorem  1  gives 

oo  oo  r* 

C  ■  /  f  P0exp[-(2ba)-1  (x8+ya)  ]  p(k/o,R/0)+Qaf8exp(-cRa/a'a) 

_00  loo  I 

.^l-p(ka/a,R(Udf) -1)  ^jdxdy  (4.1) 


where 


=  p0i3  +  (P0C/aa)  (I4-IB) 

CO 

I3  =  expt-(xa+ya)/(2ba)  ]p(k/a,R/C)dxdy 

I*  -  exp[-(xa+ya)/(2ba)-c{(x-X)a+(y-y)S}/a'a]dxdy 

IB  =  J®  exp[-(xafya)/(2b8)-cRs/aa]p(ka/a,  R(OQf)  _1)  dxdy 


(4.2) 


First, 


I4  -  exp[- (x^y8)  c/o?a]  Jg,  v,a  exp[-^(x8+y8)/(2Q'aba)+2cXx/ofa 
a 

+2cYy/o'  Jdxdy 
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=2{a3/ea)rrbsexpt-cDs/aa+2b8caD8/(aa3a)  ]  (4.3) 

using  (2.5)  of  [6],  and  where 

D8  =  Xa+Y8  as  in  Fig.  2. 

In  order  to  evaluate  I3  and  I5  we  require  a  Bessel  function,  and 
then  to  apply  a  property  of  p(R,r)  .  First,  we  rotate  the  axes  so  that 
(X, Y)  lies  on  the  new  x'-axis,  and  then  translate  the  origin  to  (X,Y)  . 
Thus  let 

x  =  (x'+DJcosoc1  -  y'sino'1 

y  =  y'coso'1  +  (x'+D)  sinor1 

where  tan  or'  =  Y/X  . 

R  is  invariant  under  translation  and  rotation,  and  Ra  =  x,a+y,a. 

So  I3  =  r  f®  exp[  (-l/(2b8) )  (x,8+y,8+2x'DfD8)  ]p(k/a,R/a)dx'dy' 

_00  '_00 

Finally,  introduce  polar  coordinates, 
x '  =  Rcos9 

y1  =  RsinG. 

Then 

I3  =  J?fZnexp[-<R*+D’)/(2b°)  -DRcose/ba]p(k/a,Va)  RlGdR 

=  Sirf exp[-(Ra+D8)/(2b8)  ]  I0  (DR/b8)  p(k/CT,R/0)  RdR,  (4.4) 

*  0 

where 

I0  (z)  =  (2tt)  1  J^8TT  exp(-zcos  0)  d  9 
is  the  modified  zero-order  Bessel  function  of  the  first  kind. 

We  require  the  result 

p(r0/s,d/s)  =  a“«r  p(rc/6,r/6)  exp[-  (r8+d8)/(2a8)  ]l0  (rd/O8)  rdr 

where  (4.5) 

s8  =  a8+68 

(See,  for  example,  Egs.  (15)  and  (16)  of  [5]). 
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Then  it  immediately  follows  from  (4.4)  and  (4.5)  that 

I3  =  an^pIk^-war^D^+C8)-15]  (4.6) 

Next,  we  make  the  same  translation  and  rotation  for  IB  as  for  I3  ,  and 
similarly  transform  to  polar  coordinates.  We  get 

IB  =  2nJ^p[ko'/C,R/(0a)  ]exp[-(Ra+D^/(2ba)-cRs/aa] 

.Io  (DR/ba)  RdR 

In  order  to  apply  (4.5)  to  this  integral,  it  is  necessary  to  rearrange 
with  the  use  of  auxiliary  factors.  Writing  IB  in  the  form 

IB  =  2nexp(-cDa/ga)J^>pJ^kaa/CJa')  ,R/(0Df)J 

.exp£-{Ra+(Dcy8/P8)  9V(2baa8/B8J.lo^D(a8  P8)  R/(bcr/P)  RdR, 

We  can  apply  (4.5)  to  get 

IB  =  2rr(Qfa/Pa)  exp(-cDa/Pa)  p[kor(aa+ba/Pa)  ^DtrP  8(aa+b8/pa) 

(4.7) 

Collecting  (4.3),  (4.6)  and  (4.7),  C  is  obtained.  We  state  this  in 
the  form  of 

Theorem  2 .  Under  the  model  described  (viz.  urban  pojwlation  density 
(4.1)  ,  damage  function  (1.1)  ,  aiming  error  as  in  Theorem  1)  ,  then  the 
expected  number  of  casualties  from  a  single  weapon  aimed  at  (X,Y)  is 
C  -  2TTP0ba£p(k/  (b8+<ja)  (ba+C3)  **)  +  (Q/p8) exp(-cDa/3a) 

.{l-Pfkap/tba+oap8)  ^Dfa/P) /(b^a3^)  S}j  (4.8) 

where  or*  *  l+2c0a ,  a>0 

p8  *  aa+2bac ,  P>o . 

Proof:  This  follows  from  (4.2),  (4.3),  (4.6)  and  (4.7). 

if  s0a  =  ba+a8,  s0>o, 

=  aa+2bac,  ^>0, 


Q.E.D. 
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Chen  (4.8)  can  be  written 

C  =  2TTP0ba|p(k/S0,D/S0)  +  QP“aexp(-cD*/38) 

.{l-pOwp/Si^/OSi))^.  (4.9) 

The  proportion  of  the  city  population  which  falls  casualty  to  this 
single  weapon  attack  is  C/(2rrPoba)  . 

The  programming  of  (4.9)  can  be  done  by  Wegner's  approximations 
[8]  to  p(R,r)  ,  and  these  are  set  out  in  the  author's  earlier  paper  [6]. 
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1  In,  a»- aarlier  paper  [A3  a  city  distributed  uniformly  over  a  circle  was  consider¬ 
ed  in  a  model  for  an  n-weapon  attack  using  the  circular  normal  damage  function  and 
aiming  error  developed  by  Hunter  [43  .  The  expected  number  of  casualties  was  evalu¬ 
ated  using  the  Circular  Coverage  Function, p(R, r)  ( [1] ,  [2] ) . 

Many  damage  curves  do  rot  fit  the  single-  or  two-term  circular  normal  form  of 
Hunter,  or  the  elliptical  normal  form  discussed  by  Grubbs  [3]  .  These  "cookie-cutter" 
curves  are  better  fitted  by  the  damage  function 

1  ,  0  <  r  <_  k 

Qexp(-cr2)  ,  r  >  k,  Q  >_  0,  c  >  0,  d-D 

where  r  is  the  distance  from  ground  zero. 

_>We  consider  a  single-weapon  attack  on  a  city  having  a  circular  normal  distribu¬ 
tion  of  population  about  the  origin  4-see  Sherratt  [7W-;  the  aiming  error^Asr; in  Hun¬ 
ter's  model  [4],  is  assumed  to  be  circular  normal  about  the  point  of  impact , 

and  the  damage  function  of  the  form  (1.1).  '  The  Circular  Coverage  Function,  diaeusaed 
in  [6],  will  again  figure  in  the  evaluation  of  the  expected  number  of  casualties. 

Section  2  evaluates  the  expected  casualties  when  the  impact  point  is  given,  and 
in  Section  3  an  expression  for  the  damage  function,  modified  to  include  aiming  error,' 
is  obtained.  Section  4  estimates  the  expected  casualties  when  the  aiming  point 
(Xi,Yi>  is  given^  To  the  writer's  knowledge,  no  previous  model  has  led  to  such  a 
result  for  the  cookie-cutter  damage  function  (1.1). 


